Helical Metal Inside a Topological Band Insulator 
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Topological defects, such as domain walls and vortices, have long fascinated physicists. A novel 
twist is added in quantum systems like the B-phase of superfluid helium He3, where vortices are 
associated with low energy excitations in the cores. Similarly, cosmic strings may be tied to propa- 
gating fermion modes. Can analogous phenomena occur in crystalline solids that host a plethora of 
topological defects? Here we show that indeed dislocation lines are associated with one dimensional 
fermionic excitations in a 'topological insulator', a novel band insulator believed to be realized in 
the bulk material Bio.gSbo.i. In contrast to fermionic excitations in a regular quantum wire, these 
modes are topologically protected like the helical edge states of the quantum spin-Hall insulator, 
and not scattered by disorder. Since dislocations are ubiquitous in real materials, these excitations 
could dominate spin and charge transport in topological insulators. Our results provide a novel 
route to creating a potentially ideal quantum wire in a bulk solid. 
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Motivated by applications to spintronics, recent theo- 
retical work on the effect of spin orbit interactions on the 
band structure of solids predicted the existence of novel 
'topological insulators' (TIs) in two[l|, Q, Q and three 
dimensions [1, [H, [f| . In two dimensional TIs (or quantum 
spin-Hall insulators), in contrast to the gap in the bulk, 
there are gapless modes at the edge which appear in time 
reversed pairs with opposite spin and velocity. Transport 
measurements on HgCdTe quantum wells [ll H| have pro- 
vided evidence for the existence of these 'helical' edge 
states[H,i,[i3. In three dimensions, topological insula- 
tors are classified as strong (z/q = 1) or weak (vq = 0), the 
former has surface states with an odd number of Dirac 
points. In addition, 3D TIs are characterized by indices 
(y\, V2, v%) with respect to a basis of reciprocal lattice 
vectors {G\, G 2 , G3). They define a time reversal invari- 
ant momentum (TRIM) M v = \[y x G x + v 2 G 2 + ^3<5 3 ). 
The recent prediction Q and evidence from angle re- 
solved photoemission experiments [H| for a strong TI 
phase in the alloy Bio.gSbo.i, has led to heightened inter- 
est in these systems. 




FIG. 1: (color online) [a] Diamond lattice screw dislocation 
along B — a(l, 1,0) direction (shown as a cylinder with red 
ends). The two diamond sublattices are shown as yellow and 
white spheres. A cutting plane P that was used to generate 
the dislocation is shown, orthogonal to the strength t bond 
along o/2(l,— 1,-1). [b] Screw dislocation in a stacked 2D 
topological insulator. Unpaired edge modes (one member of 
the edge mode pair is shown in red) of the top and bottom 
layer must propagate through the dislocation by continuity. 



Dislocations are line defects of the three dimensional 
crystalline order, characterized by a vector B (Burgers 
vector) , which is a lattice vector. It must remain constant 
over the entire dislocation. The center of the dislocation 
is at R(cr) where a parameterizes the line defect. A con- 
venient way to visualize a dislocation is via the Volterra 
process. One begins with the perfect crystal and chooses 
a plane P that terminates along the curve R(a) where 
the defect is to be produced. The crystal on one side of 
the plane P is then displaced by the lattice vector B, and 
additional atoms are inserted or removed if required. At 
the end, crystalline order is restored everywhere except 
near the curve R(a). A screw dislocation (see Figured]) is 
a line defect with its tangent vector t = dR(cr)/da \\ B , 
while an edge dislocation is also a line defect with t _L B. 
In general, a dislocation varies between these two simple 
types along its length. 

We now state the main result of this paper. Consider a 
dislocation with Burgers vector B inside a three dimen- 
sional topological insulator characterized by (vq; M„). Iff 



B ■ M v = 7r (mod 2tt) 



(1) 



the dislocation induces a pair of one dimensional modes 
bound to it that traverse the bulk band gap. The modes 
are related by the Kramers time reversal symmetry trans- 
formation and reflect the nontrivial topology of the band 
structure. We shall refer to this one dimensional state 
as the 'helical metal'. While a trivial insulator could 
also develop one dimensional propagating modes along 
a dislocation, the precise count obtained here - a single 
Kramers pair - cannot arise. The topological stability of 
the helical metal results from time reversal symmetry, 
which prohibits backscattering between the oppositely 
propagating modes 0. From Eqn. [T]we note that some 
dislocations e.g. with B _L M v do not carry these heli- 
cal modes. Also, there is the 'pristine' strong topological 
insulator vq — 1, M v — 0, where none of the disloca- 
tions induce helical modes. Experiments indicate that 
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Bio.gSbo.i is a strong TI with a nontrivial M u . 

Dislocations in the Diamond Lattice TI Instead 
of plunging into a general proof of our results, we first 
present a numerical calculation of a dislocation in a sim- 
ple TI model on the diamond lattice 

8A 

H = *I>*W + *7^f E 4(4 x 4) • <wc feCT ,(2) 

where d^, are the two nearest neighbor bond vec- 
tors leading from site i to fc, er are the spin Pauli ma- 
trices, and 2a is the cubic cell size. The system is a 
Dirac metal at this point, and a full gap is opened upon 
distorting the lattice, by making the nearest neighbor 
hopping strengths t + St along one of the directions e.g. 
|(1,1,1). If St > (St < 0), the system is a strong 
(weak) topological insulator with v = 1 (y = 0) and 
M v = £(1,1,1). We now study a screw dislocation 
and use periodic boundary conditions - which requires 
us to consider a separated pair of dislocations. Transla- 
tion symmetry is present along the dislocation axis, and 
we label states by k, the crystal momentum in this di- 
rection. The precise geometry we use is shown in Figure 
QJi. If a*i = a(0, 1, 1), a 2 = a(l, 0, 1), — a(l, 1, 0) are 
the basis vectors, the dislocations are along the 03 axis. 

For a pair of unit screw dislocations with B = a(l, 1, 0), 
we find four one dimensional modes, two down moving 
and two up moving modes, within the bulk gap (Figure 
[5^). The wavefunctions of these modes are peaked along 
the two dislocations (Figure [2}:) - hence a pair of oppo- 
sitely propagating modes is present for each dislocation 
line. The modes are time reversal (Kramers) conjugates 
of each other. Not all dislocations carry these gapless 
modes - if we double the Burgers vector of this disloca- 
tion (Figure^), or consider a different relative orienta- 
tion when the unequal bond is orthogonal to the Burgers 
vector (Figure [2ji), the one dimensional modes do not 
appear. The sign of St is immaterial - both weak and 
strong topological insulators show this physics. 

Note the condition for a gapless line mode to exist on 
a dislocation cannot depend on its orientation, since the 
modes must propagate through the entire defect. Hence, 
it can only depend on the Burgers vector B. Therefore, 
although we have explicitly considered the case of a screw 
dislocation, the result also holds for an edge or mixed 
dislocation with the same Burgers vector. 

We now present an analytical rationalization for the 
existence of gapless modes along the dislocation line, in 
this model of the strong TI. We study the model in a limit 
that is analytically tractable, but nevertheless adiabati- 
cally connected to the parameter regime of interest. Con- 
sider the imaginary plane P (orthogonal to (1,-1,-1)), 
used to construct the dislocation, which ends along the 
dislocation line. Extend this plane through the entire 
crystal, weakening all bonds that are cut by it. Now, the 
crystal is almost disconnected into two disjoint pieces. 




FIG. 2: (color online) Electronic spectrum of the diamond 
lattice strong TI (M„ = £(1, 1, 1)) in the presence of a pair 
of screw dislocations. Electronic states are shown as a func- 
tion of k, the wavevector along the dislocation. For [a][b][c], 
the dislocations are oriented along the S3 — o(l, 1, 0) axis, [a] 
Spectrum, when the dislocation Burgers vector B = S3 and 
satisfies B-M v — it. A pair of counter-propagating modes per 
dislocation is found, that span the band gap (each midgap 
level happens to be doubly degenerate), [b] Spectrum when 
B — 2S3 and B ■ M v = Omod 2ir. No dislocation modes occur, 
[c] Probability distribution of a pair of midgap modes shown 
in the Si, S2 plane that intersects the dislocations. The modes 
are bound near the dislocations, [d] The dislocations are ori- 
ented along B = a(0, 1, —1) and B ■ M v — 0. No dislocation 
modes occur. Calculations were done on a 36x36x18 unit cell 
system with periodic boundary condition along Si, S2, B di- 
rections, and dislocations separated by a half system size. The 
hopping parameters used in Eqn[2]are t = 1, t + 5t = 2, and 
Xso = 0.125. 

Gapless states are generated on the two surfaces S+, S-, 
each with a single Dirac node located at the surface mo- 
mentum rfiD — j-(l, 1,0). Electronic states in the vicin- 
ity of the Dirac point can be expanded as i b{xi,X2) ~ 
e lrn,D ' x if}(x), where i/j(x) varies slowly over a lattice spac- 
ing. The effective Hamiltonian for such states on the 
two surfaces is ±Hq = p±a x + P2& y [1, [13], where pi, 
P2 are surface momenta measured from the Dirac point, 
with pi along the dislocation line (the (1,1,0) direction). 
Thus if Hz = ±1 represents the two surfaces, the effec- 
tive Hamiltonian is H = H+(p)/i z . Now, if we imag- 
ine reconnecting the surfaces without creating a disloca- 
tion, we add a hopping term between the two surfaces: 
H = H^[p)fx z + mfJ, x , which leads to fully gapped disper- 
sion, with gap 2m, which can be adiabatically connected 
to the uniform bulk insulator. 

Now consider introducing a screw dislocation, by dis- 
placing the surface to one side of the plane P by the 
lattice vector B = o(l,l, 0). In our surface coordi- 
nate system, the dislocation line is along the x\ axis, 
at X2 — 0, and the plane P is at X2 > 0. The crucial 
observation is that when reconnecting the bonds across 
the plane P, if the Dirac node is at momentum mjj with 
mjj ■ B = 7T, then the Dirac particles acquire a phase 
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shift of 7r. The effective Hamiltonian in the presence of 
a dislocation is then: Ha ls = H + (—ih\7)n z + m(x2)\i x 
where the Dirac mass term is: m(x 2 < 0) = m but 
m(x2 > 0) = — m. Such Dirac Hamiltonians are well 
known to lead to low energy modes [13j that are protected 
by index theorems. At pi = 0, a pair of midgap modes 

are present ip±(x) = e 7 ^ J° d3/m ^^>o,± where ipo,± are 
the two solutions of: fi y a y ipo = i/jq. For finite p\ these 
split to give rise to the left and right moving modes. 

More generally, if there are several such surface Dirac 
modes m % D , what is required is that an odd number of 
them acquire a n phase shift on circling the dislocation 
Y^i ^D ' B = 7r(mod27r). This guarantees the existence 
of at least a pair of gapless ID modes on the dislocation. 
If on the other hand, ^ rn l D ■ B = 0(mod 2tt), there are 
no protected modes on the dislocation. 

Dislocations in a General TI Armed with these in- 
sights for a specific model, one may discuss the general 
case of an arbitrary dislocation in a topological band in- 
sulator. The condition for the existence of a protected 1 
D mode is given by equation [TJ and derived in the gen- 
eral case in the paragraph below. For the weak TI case, 
however, a more intuitive derivation of this result is pos- 
sible. The weak TI is adiabatically connected to a stack 
of decoupled two dimensional TIs, stacked along the M v 
direction. The top and bottom surfaces of the stack are 
fully gapped. Consider creating a screw dislocation in 
such a stack, by cutting and regluing layers. Cutting a 
plane results in a pair of edge modes, which are gapped 
after the planes are glued back together. However, edge 
modes in the top and bottom layers are left out in this 
process (see Figure [Tp). Since a single pair of such modes 
cannot begin or end, the helical mode must propagate 
through the dislocation core to complete the circuit. A 
similar result holds for three dimensional Chern insula- 
tors [1J|, characterized by a reciprocal lattice vector G . 



TRIMs mi, TO2 is a topologically protected quantity jij: 



There, a dislocation line should have G °' B chiral modes 
propagating along it. 

We now derive the main result equation [TJ First, let 
us briefly review the meaning of the topological insu- 
lator invariants (uq; M v ). These are conveniently ex- 
pressed in terms of quantities £>i = („ in2 „ 3 ) with rii = 0, 1 
defined at the 8 TRIMs of the 3D Brillouin zone(BZ) 

- i=(n 1 n 2 n 3 ) 



(niGi + 712G2 + H3G3Y/2, which can be 



computed for a given band structure jj]. The gauge in- 
variant indices are given by: (— l) Ua = Y[ n = q 1 ^n 1 n 2 n 3 



and (-l)"* 



n, 



The indices 



Vi carry information about the surface band structure, in 
particular, the number of times that the surface bands 
cross a generic Fermi energy which lies within the bulk 
gap, along a line in the surface BZ. For example con- 
sider the surface spanned by the basis vectors Si , a 2 nor- 
mal to G3. The parity (even vs odd) of the number of 
such band crossings N cross , when connecting two surface 



(-1) 



AY 



1 ^rfi.o ^ - 

mi- 



(3) 



Thus v\ = iV cross (mod 2) between fh\ = and % 



Gl + G2 on this surface. 

Without loss of generality, consider now creating a 
screw dislocation in a TI with Burgers vector B = a\. 
This is created using the Volterra procedure by a cutting 
3.1 — 0,2 surface containing the dislocation. Atoms on one 
side of the surface are displaced by B and reconnected. 
The crucial step in our analytical derivation of the helical 
modes was whether there are even or odd number of sur- 
face Dirac nodes that acquire a 7r phase shift on creating 
the dislocation. Although in general, the surface band 
structure may be quite complicated, for the purpose of 
deriving robust topological properties one may adiabati- 
cally deform the crystal structure (during which the bulk 
gap remains open, and the Bravais lattice structure stays 
fixed) so that the surface modes are always Dirac nodes 
that cross the band gap, centered at the surface TRIMs. 

}for 



Then, there are two surface TRIMs : 

which B ■ m = 7r(mod27r), i.e. acquire a 7r phase shift on 
creating the dislocation. If the number of surface Dirac 
nodes present in total at these two TRIMs A^Dirac is odd, 
then the dislocation will host a protected helical mode. 
Clearly, the parity of iVnirac is the same as the number 
of band crossings between the two surface TRIMs given 
by Equation[3J Thus, the condition for the helical modes 
is vi = 1 or M u ■ B — 7r(mod27r). 

In general, the Burgers vector is a multiple of the prim- 
itive lattice vector. If it an odd multiple, the derivation 
above is unaffected. If it is an even multiple, there are no 
dislocation modes since the phase shifts are always triv- 
ial. Hence we arrive at [TJ for the existence of protected 
helical modes along the dislocation. 

Effect of Disorder It is sometimes stated that in the 
presence of disorder, only the Uq index is robust, while 
the M„ index is irrelevant. Hence weak TIs are believed 
to be unstable in the presence of disorder 0] ■ Since our 
results are controlled precisely by the M v index, they 
shed light on when they retain meaning. If the disorder 
is so strong that the dislocations are no longer well de- 
fined, then the only remaining distinction is the vq index. 
Similarly, an insulator with M v ^ can be converted to 
one with M v — 0, by introducing a potential modulation 
with period M v . Elementary dislocations that carry he- 
lical modes will then cost infinite energy per unit length. 
However, for a sufficiently weak disorder potential, the 
dislocations arc expected to remain well defined, in which 
case the distinction captured by M v is physically rele- 
vant. Certainly, in real crystals which are inevitably dis- 
ordered, dislocations are well defined objects, and hence 
the physics discussed in this paper should apply. 

Experimental Consequences In real solids, disloca- 
tions are always present and the predicted helical modes 
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should have important experimental consequences. Scan- 
ning tunneling microscopy (STM) of a topological insu- 
lator surface where dislocations terminate, is particularly 
well suited to verify these predictions. Since the precise 
atomic arrangement is visualized by STM, the nature of 
the dislocations involved can be characterized. At the 
same time, the finite density of states associated with 
the one dimensional modes will lead to an enhanced tun- 
neling density of states near the dislocation core. This 
will be particularly striking if one tunes to an energy 
where the surface modes contribute minimally (eg. at 
the Dirac node). Such experiments are immediately fea- 
sible on the putative topological insulator Bio.gSno.i with 
A7 structure. Since the band structure of that material 
is complex and still debated [lq . we demonstrate quali- 
tatively what is expected, in the simpler diamond lattice 
strong TI model [2] Moreover, since the crystal symme- 
try in the two cases share several features - such as a 
three fold axis parallel to the M„ vector, some predictions 
for the diamond lattice model may be directly relevant 
to the A7 structure. First consider the surface orthog- 
onal to the strong bond (1,1,1) direction. We choose 
this surface to cut three strength t bonds (rather than 
a single strong bond)in order to obtain a surface band 
structure similar in some respects to the (111) surface 
of Bio.gSno.i. We call this the (HI)' surface. This sur- 
face has a single Dirac node centered at the T point. 
We consider a pair of separated screw dislocations with 
B = a(l, 1,0), which carry ID helical modes, that ter- 
minate on the surface. STM measures the local density 
of states (LDOS) at a given energy and surface location 
p(f,E) = J2 a \' l l } a{r)\ 2 5{E — E a ), where the a sum runs 
over all eigenstates. The dislocations appear as peaks 
in the LDOS as shown in Figure [3^,, due to these heli- 
cal modes. Note, the the ledge connecting the two screw 
dislocations on the surface also shows an enhanced DOS. 

More striking evidence for the relation [T] appears when 
we consider the edge dislocations on the (111) surface. 
If the Burgers vector of the dislocation is B — a{\, 0, 1) 
the one dimensional modes are present and visible in the 
LDOS (Figure[3b), but if B = a(l, -1,0) then they are 
absent (Figure . Here, there is no physical line con- 
necting the two defects. 

Since elastic scattering by nonmagnetic impurities can- 
not scatter electrons between the counter-propagating 
helical modes of the dislocation, they behave as ideal 
quantum wires which may be relevant for spin and charge 
transport applications at low temperatures when inelas- 
tic processes are frozen out. The challenge will be to 
separate this conduction mechanism from protected sur- 
face mode conduction, which is also present. In conduc- 
tion across the short direction of anisotropic samples (eg. 
disc shaped), dislocation induced direct conduction paths 
should dominate over surface conduction. If dislocations 
can be induced in a controlled fashion during the growth 
process, then the direction dependent condition for the 
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FIG. 3: (color online) Surface LDOS as measured by STM, 
for various dislocation configurations, in the diamond lattice 
strong TI with parameters as in Fig. 2. In all cases disloca- 
tion pairs are separated by a half system size, and are di- 
rected along S3. Topography of the integrated LDOS in the 
energy window [—0.1, 0.1] is shown, with insets displaying en- 
ergy dependent LDOS at special points on the surface, [a] 
The (111)' surface intersecting a pair of screw dislocations 
with B = o(l,l,0). These carry ID modes which give rise 
to a peak in the integrated LDOS. In [b] and [c], edge dis- 
locations on the (111) surface. In [b], the dislocation sat- 
isfies B ■ M v — tv, and has ID modes visible in the LDOS 
while in [c] the dislocation has B ■ M„ = 0, hence no ID line 
mode and no enhanced LDOS. In [a][b][c] only the top layer of 
atoms is shown, and calculations are done on 36x36x18 unit 
cell system, with periodic boundary conditions in the 01, S2 
directions. 



existence of helical modes, and dislocation density de- 
pendence of conductivity can be used to isolate this con- 
tribution. Another possibility is to introduce magnetic 
impurities on the surface of the insulator, which could 
potentially localize surface modes but have little impact 
on dislocation helical modes deep in the bulk. A rough 
estimate of the excess electric conductivity induced by 
dislocation modes for a dislocation density ~ 10 12 m~ 2 
and a low temperature scattering length / ~ lpm along 
the dislocation [8| yields p = ~ lOmfhii, which 

should be experimentally accessible. 
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